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ABSTRACT 

We explore the quenching of low-mass galaxies (10"^ < M^ ^10^ Mq) as a function of lookback time using 
the star formation histories (SFHs) of 38 Local Group dwarf galaxies. The SFHs were derived from analyz¬ 
ing color-magnitude diagrams of resolved stellar populations in archival Hubble Space TelescopeAVide Field 
Planetary Camera 2 imaging. We find: (1) Lower mass galaxies quench earlier than higher mass galaxies; (2) 
Inside of Rymai there is no correlation between a satellite’s current proximity to a massive host and its quenching 
epoch; (3) There are hints of systematic differences in quenching times of M31 and Milky Way (MW) satellites, 
although the sample sample size and uncertainties in the SFHs of M31 dwarfs prohibit definitive conclusions. 
Combined with literature results, we qualitatively consider the redshift evolution (z = 0 - 1) of the quenched 
galaxy fraction over ^ 7 dex in stellar mass (lO"^ < M^ < Mq). The quenched fraction of all galaxies 

generally increases toward the present, with both the lowest and highest mass systems exhibiting the largest 
quenched fractions at all redshifts. In contrast, galaxies between M^ ^ 10^-10^^ M© have the lowest quenched 
fractions. We suggest that such intermediate-mass galaxies are the least efficient at quenching. Finally, we 
compare our quenching times with predictions for infall times of low-mass galaxies associated with the MW. 

We find that some of the lowest-mass satellites (e.g., CVn H, Leo IV) may have been quenched before infall 
while higher mass satellites (e.g., Leo I, Fornax) typically quench ^1-4 Gyr after infall. 

Keywords: galaxies: dwarf - galaxies: evolution - galaxies: stellar content - Hertzsprung - Russell and C - 
M diagrams - Local Group 


1. INTRODUCTION 

Environment is believed to play an important role in the 
evolution of low-mass galaxies. Their shallow potentials 
make them vulnerable to a variety of internal (e.g., galac¬ 
tic winds, supernovae feedback) and external (e.g., tidal ef¬ 
fects, ram pressure stripping) processes that can reshape and 
remove their baryons. The central role of environment is ob- 
servationally borne out by the dwarf ‘morphology-density’ re¬ 
lationship, which shows that red, gas-poor, non-star forming 
dwarf galaxies are predominantly found in close proximity 
to a massive host, while blue, gas-rich, star-forming dwarfs 
are preferentially located in lower density environments (e.g., 
Hodge 1971; Einasto et al. 1974; van den Bergh 1994; Mateo 
1998; van den Bergh 2000; Grcevich & Putman 2009; Tol¬ 
stoy et al. 2009; Weisz et al. 2011; Geha et al. 2012; Mc- 
Connachie 2012; Spekkens et al. 2014). This powerful obser¬ 
vational constraint favors models in which quenching (i.e., the 
complete cessation of star formation)^ is induced by environ- 
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mental mechanisms (e.g., Mayer et al. 2001,2006) as opposed 
to purely internal processes (e.g., Dekel & Silk 1986). 

Despite the clear importance of environment for gas re¬ 
moval and quenching in low-mass galaxies, our knowledge 
of the responsible mechanisms is less certain. Tidal effects 
and ram pressure are believed to be the primary mechanisms 
for shutting down star formation (e.g., Gunn & Gott 1972; 
Einasto et al. 1974; Faber & Lin 1983^ However, their rel¬ 
ative effectiveness is modulated by a number of factors in¬ 
cluding halo mass and accretion history, infall time, proxim¬ 
ity to a host, stellar feedback (which can induce outflows and 
weaken the effective potential), and the heating of gas from 
cosmic reionization in the lowest mass systems (e.g., Dekel 
& Silk 1986; Thoul & Weinberg 1996; Gnedin 2000; Mayer 
et al. 2001; Busha et al. 2010; Govemato et al. 2010; Zolotov 
et al. 2012; Lu et al. 2014; Wetzel et al. 2015). Because the 
complex interplay between each of these processes can vary 
with time, the morphology-density relationship at z = 0 cannot 
strongly discriminate between various scenarios for creating 
environmentally dependent quenching in dwarf galaxies. 

Orthogonal constraints on quenching mechanisms come 
from star formation histories (SFHs). By tracing star for¬ 
mation over cosmic time, it is possible to determine exactly 
when and how quickly a dwarf galaxy quenched, providing 
insights into the physical mechanisms responsible for shut¬ 
ting down star formation. SFHs for galaxies located within a 
few Mpc can be reconstructed in detail by analyzing their re¬ 
in these systems, the re-ignition of star formation is unlikely without signifi¬ 
cant gas accretion. We view such quenched systems as qualitatively different 
from post-starburst galaxies, despite the use of similar ‘quenching’ terminol¬ 
ogy in the literature. Generally, post-starburst galaxies have had recent but 
no current star formation (i.e., ultra-violet flux but no detectable Ha), and 
typically have large gas reservoirs, suggesting that they may be capable of 
restarting star formation in the future. 
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solved star color-magnitude diagrams (CMDs; e.g., Tosi et al. 
1989; Harris & Zaritsky 2001; Dolphin 2002; Aparicio & Hi¬ 
dalgo 2009). The power of SFHs for constraining quenching 
was recently illustrated by Sohn et al. (2013) in their analysis 
of Milky Way (MW) satellite Leo I. The reconstructed orbital 
history has revealed that Leo I entered the virial radius of the 
MW 2 Gyr ago and made a pericentric passage ^ 1 Gyr ago 
at a distance of 90 kpc. The SFH shows an enhancement fol¬ 
lowed by complete cessation of star formation over the same 
timescale, suggestive of a rapid quenching scenario likely due 
to ram pressure stripping, as its pericentric distance seems too 
large for tidal effects to be significant. 

Despite the clear utility of SFHs for quenching studies, 
there have been few efforts to systematically characterize and 
analyze quenching in nearby dwarfs. The majority of rele¬ 
vant studies largely focus on quantifying the effects of early 
reionization on the lowest mass (M^ < 10® M 0 ) MW satel- 
lites (e.g.. Sand et al. 2009,2010; Okamoto et al. 2012; Brown 
et al. 2012,2014). Outside of the MW sub-group, there is only 
a small sample of quenched, low-mass galaxies with well- 
constrained SFHs (Cetus, Tucana, Andromeda II, Andromeda 
XVI; Monelli et al. 2010a,b; Weisz et al. 2014c), making it 
challenging to explore broad quenching trends in low-mass 
systems. At the same time several studies have compared 
N-body simulations to z = 0 properties of low-mass galaxies 
to constrain quenching timescales (e.g.. Slater & Bell 2013, 
2014; Phillips et al. 2014; Wheeler et al. 2014; Phillips et al. 
2015), but combining SFHs with the simulations is still an 
emerging field (e.g., Rocha et al. 2012). 

In this paper, we use the SFHs of 38 LG dwarf galax¬ 
ies to systematically measure the quenching epochs and 
timescales in low-mass systems (M^ ^10^ Mq). The SFHs 
are taken from the Weisz et al. (2014a) uniform analysis 
of optical CMDs based on archival imaging taken with the 
Hubble Space Telescope / Wide Field Planetary Camera 2 
(HSTAVFPC2; Holtzman et al. 1995). Using these SFHs, 
we investigate empirical trends in quenching epochs and 
timescales as functions of basic galaxy properties (e.g., stel¬ 
lar mass, proximity to a host), place the LG into the context 
of broader studies of quenching, and compare our data with 
select models of satellite infall in the LG. We have previously 
used this dataset to examine the observational evidence for 
reionization as a quenching mechanism in the lowest mass 
systems, and refer readers to Weisz et al. (2014b) for a more 
detailed discussion of this specific quenching mechanism. 

This paper is organized as follows. We summarize the sam¬ 
ple selection, photometry, and SFH measurement technique 
in §2. In §3, we explore empirical trends in the fraction of 
quenched galaxies as a function of mass and environment, 
and place our results in the context of the broader universe 
using data from the literature. In §4, we quantify delay times 
between infall and quenching by comparing our results with 
the model analysis presented in Rocha et al. (2012). Finally, 
we summarize our findings in §5. Throughout this paper, the 
conversion between age and redshift assumes the cosmology 
as detailed in Planck Collaboration et al. (2014). We recog¬ 
nize that several definitions of quenching epoch can be found 
in the literature. For direct comparison, we have provided all 
SFH data in tabulated digital form in Weisz et al. (2014a). 


2. THE DATA 

2.1. Galaxy Sample and Star Formation Histories 


Our sample contains 38 LG dwarf galaxies that all have 
deep archival HSTAVFPC2 imaging. The sample includes 
25 diverse quenched galaxies (dwarf spheriodals, dSphs; and 
dwarf ellipticals, dEs) spanning a large dynamic range in stel¬ 
lar mass (10^ < M^ ^10^ Mq). Our sample also contains 
5 transition dwarfs (dTrans), which have HI but no evidence 
of on-going star formation (i.e., no detectable Ha; e.g., Ma¬ 
teo 1998), and 8 star-forming dwarfs (dwarf irregulars, dirrs), 
both of which provide good control and comparison samples. 
The properties of our galaxy sample are listed in Table 1 . 

We selected galaxies based on the availability of deep 
archival HSTAVFPC2 observations. All photometry was car¬ 
ried out using HSTPHOT^, a photometry package designed 
specifically for the under-sampled point spread function of 
WFPC2 (Dolphin 2000), as part of the Local Group Stellar 
Populations Archive project^^ (LOGPHOT; Holtzman et al. 
2006). To characterize completeness and observational biases 
in the photometry, we ran ^10^ artificial star tests per HST 
field, i.e., we inserted stars of a known magnitude into our 
images and recovered their photometry in an identical man¬ 
ner to the photometry of real stars. Details of the photometric 
reduction process can be found in Holtzman et al. (2006). 

All SFHs used in this paper were derived using the CMD 
fitting package MATCH (Dolphin 2002), as discussed exten¬ 
sively in Weisz et al. (2014a); we refer the reader to that pa¬ 
per for full details. Throughout the paper, we quote the 68% 
confidence interval of the total uncertainties, including both 
random (i.e., statistical uncertainties) and systematic (i.e., a 
proxy for variations in the SFH as if they were solved with 
different stellar models) components. Details of the uncer¬ 
tainty calculations for this sample are presented in Weisz et al. 
(2014a), and the general methodological details of the sys¬ 
tematic and random uncertainty determinations are presented 
in Dolphin (2012) and Dolphin (2013), respectively. 

2.2. Defining an Epoch of Quenching 

As a first step, we must define a quenching epoch for galax¬ 
ies in our sample. Identifying the quenching epoch for such 
a diverse collection of galaxies using CMD-based SFHs in¬ 
volves some subtlety, primarily owing to the presence of blue 
straggler stars (BSs). BSs — believed to be the result of 
merged or extreme mass transfer among low-mass main se¬ 
quence stars (e.g., Davies 2015, and references therein) — 
have similar optical colors and luminosities as intermediate 
age main sequence stars (2-5 Gyr old). Therefore, they may 
be mistaken as a signpost of intermediate age star formation, 
particularly in lower mass, quenched galaxies, where the ‘blue 
plume’ is typically sparse. 

Unfortunately, accounting for BSs in CMD-based SFH 
measurements is non-trivial. Because there are no generative 
models for BSs, it is impossible to explicitly include them 
in the CMD fitting process. One often-used solution is to 
simply exclude putative BSs from a CMD when measuring 
a SFH. While this may be suitable for galaxies in which there 
is a strong prior belief that the ‘blue plume’ is composed of 
only BSs (e.g., the faintest MW satellites; Brown et al. 2012, 
2014). It is challenging to implement this approach in a large, 
diverse sample of galaxies. Several studies have shown that 
increasingly luminous galaxies have blue plumes that contain 
a larger fraction of genuine main sequence stars (e.g., Mapelli 
et al. 2007; Momany et al. 2007; Mapelli et al. 2009). Thus, 

^ http://purcell.as.arizona.edu/hstphot/ 
http://astronomy.nmsu.edu/logphot 
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Table 1 

Observational Properties of our Local Group Dwarf Sample 


Galaxy 

Morphological 

Mv 

M* 

Mhi 

DhosI 

fh 

T90 

Name 

Type 


(10" Mq) 

(lO** Mq) 

(kpc) 

(pc) 

log(yr) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Andromeda I 

dSph 

-11.7 

3.9 

0.0 

58 

672 

Q o^+0-0>015 
^•°Wo.05,0.14 

Andromeda II 

dSph 

-12.4 

7.6 

0.0 

184 

1176 

Q ^r+0.01,0.07 
^•'^-0.01,0.05 

Andromeda III 

dSph 

-10.0 

0.83 

0.0 

75 

479 

Q ^.^+0.05,0.3 
-0.0,0.06 

Andromeda V 

dSph 

-9.1 

0.39 

0.0 

no 

315 

Q g.,+0.02,0.23 
'-0.06,0.31 

Andromeda VI 

dSph 

-11.3 

2.8 

0.0 

269 

524 

Q .70+0.06,0.3 
^•''’-0.02,0.09 

Andromeda XI 

dSph 

-6.9 

0.049 

0.0 

104 

157 

1^11 +0.0,0.0 
^^•^^-0.08A2 

Andromeda XII 

dSph 

-6.4 

0.031 

0.0 

133 

304 

Q r^+0.2,0.33 
^•^”-0.04,0.04 

Andromeda XIII 

dSph 

-6.7 

0.041 

0.0 

180 

207 

+0.0,0.0 
iC'.r_0,44 0.46 

Carina 

dSph 

-9.1 

0.38 

0.0 

107 

250 

Q 4<+0 0A11 
^•^”-0.1,0.12 

Canes Venatici I 

dSph 

-8.6 

0.23 

0.0 

218 

564 

Q +0.0,0.06 

^•^^-0.1,0.11 

Canes Venatici II 

dSph 

-4.9 

0.0079 

0.0 

161 

74 

Q q/^+0.0,0.05 
^•^^-0.24,0.24 

DDO 210 

dTrans 

-10.6 

1.6 

4.1 

1066 

458 

Q co+0.07,0.52 
^•^^-0.16,0.37 

Draco 

dSph 

-8.8 

0.29 

0.0 

76 

221 

1 n 01 ,0.06 

iU.Ul_o 03,0.11 

Fornax 

dSph 

-13.4 

20.0 

0.0 

149 

710 

Q 0 0+0.02,0.14 
^•■^^-0.02,0.06 

Hercules 

dSph 

-6.6 

0.037 

0.0 

126 

330 

^ ^ 1 +0.0,0.0 
^'^•^-0.61,0.61 

ICIO 

dlrr 

-15.0 

86.0 

50.0 

252 

612 

^•^^-0.0,0.06 

Q q+0.0,0.08 

^•■^-0.07,0.12 

IC 1613 

dlrr 

-15.2 

100.0 

65.0 

520 

1496 

Leo A 

dlrr 

-12.1 

6.0 

11.0 

803 

499 

0 ^0+0.0,0.22 
-0.0,0.08 

Leo I 

dSph 

-12.0 

5.5 

0.0 

258 

251 

Q ^0+0 0,0.05 
^•^■^-0.0,0.0 

Leo II 

dSph 

-9.8 

0.74 

0.0 

236 

176 

Q Q1 +0.0,0.05 
^•^^-0.0,0.04 

Leo IV 

dSph 

-5.8 

0.019 

0.0 

155 

206 

^^•^'^-0.57,0.57 

Leo T 

dTrans 

-8.0 

0.14 

0.28 

422 

120 

Q .70+0.0,0.0 

^•^'’-0.02,0.05 

LGS3 

dTrans 

-10.1 

0.96 

0.38 

269 

470 

Q cq+0 03,0.19 
^•^°-0.02,0.09 

M32 

dE 

-16.4 

320.0 

0.0 

23 

no 

Q .70+0.0,0.49 
^•^■^-0.0,0.0 

NGC 147 

dE 

-14.6 

62.0 

0.0 

142 

623 

Q ^0+0.02,0.53 
^•‘^■^-0.0,0.01 

9 56+0 02,0.4 

^•'^”-0.02,0.05 

NGC 185 

dE 

-14.8 

68.0 

0.11 

187 

42 

NGC 205 

dE 

-16.5 

330.0 

0.4 

42 

590 

Q 04+0.0,0.47 
^•■^^-0.02,0.04 

NGC 6822 

dlrr 

-15.2 

100.0 

130.0 

452 

354 

Q ^+0.04,0.58 
^•^-0.0,0.0 

PegDIG 

dTrans 

-12.2 

6.61 

5.9 

474 

562 

Q 00+0.0,0.82 
-0.11,0.11 

Phoenix 

dTrans 

-9.9 

0.77 

0.12 

415 

454 

Q ^0+0.02,0.08 
^•‘^'’-0.0,0.01 

Q Q1 +0.0,0.35 
°-^^-0.06,0.3 

SagDIG 

dlrr 

-11.5 

3.5 

8.8 

1059 

282 

Sagittarius 

dSph 

-13.5 

21.0 

0.0 

18 

2587 

Q co+0.0,0.19 
^•^'’-0.04,0.04 

Sculptor 

dSph 

-11.1 

2.3 

0.0 

89 

283 

1 n ao+O-02,0.05 
-0.01,0.17 

Sex A 

dlrr 

-14.3 

44.0 

77.0 

1435 

1029 

0 qc+0.0,0.36 
°-°^-0.0,0.07 

Q 0^+0.02,0.79 
^•^”-0.03,0.08 

Q 01 +0.0,0.2 

^•°^-0.0,0.19 

Sex B 

dlrr 

-14.4 

52.0 

51.0 

1430 

440 

Tucana 

dSph 

-9.5 

0.56 

0.0 

882 

284 

Ursa Minor 

dSph 

-8.8 

0.29 

0.0 

78 

181 

-0.1,0.18 

WLM 

dirr 

-14.2 

43.0 

61.0 

836 

2111 

Q qc+0.06,0.27 
-0.0,0.0 


Note. — Physical properties for galaxies in our sample taken from McConnachie (2012). The stellar mass 
listed in column (4) is computed from the integrated V-band luminosity and assumes Mq/Lq= 1. Column (8) 
indicates the epoch at which 90% of the stellar mass formed, including the random uncertainties (first error 
term) and the total uncertainty (second error term). Both reflect their respective 68% confidence intervals. 
Note that for completeness we have also included these values for galaxies that are not currently quenched, 
i.e., dIrr/dTrans morphological types. 


excluding blue stars in all early-type systems a priori would 
mute the signatures of true intermediate age star-formation. 
Alternatively, excluding BSs in some galaxies and not in oth¬ 
ers is equally challenging, due to ambiguities in the point at 
which the ‘blue plume’ is no longer purely BSs. Further, such 
an approach would compromise the homogeneity of our anal¬ 
ysis. Thus, as discussed in Weisz et al. (2014a), we have 
elected to not account for blue stragglers in the CMD mod¬ 
eling process, and instead rely on judicious interpretation of 
the resulting SFHs. 

For the purpose of this paper, we use the lookback time at 
which a galaxy formed 90% of its total stellar mass (hereafter 
Tqo) as a proxy for the epoch at which a galaxy quenched. 
This value was selected to minimize the impact of BSs on de¬ 
termining the quenching epochs of the least luminous systems 


in our sample. The presence of BSs tends to shift the youngest 
bona fide star formation to younger ages, while keeping about 
the same total mass formed in the system. To estimate the im¬ 
pact of BSs on the SFHs, we re-ran SFHs of several lower 
luminosity dSphs having masked out the ‘blue plume’ and 
found a typical mass contribution from BSs of ~ 5%, with 
an upper limit of < 10%. Thus, while our adoption of rqo 
is a conservative choice, it does avoid potential confusion of 
the quenching epoch due to the presence of BSs. Readers who 
wish to explore other quenching metrics can electronically ob¬ 
tain all tabulated SFH data from Weisz et al. (2014a). 

We list the values of rqo for each galaxy in Table 1 . These 
values were computed by interpolating the cumulative SFHs 
in linear time for each galaxy, including both random and sys¬ 
tematic uncertainties, onto a fine time grid of 10 Myr in spac- 
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DmW,M31 (kpc) 


Figure 1. Distance to the nearest host plotted versus present day stellar mass. 
Galaxies in our sample are color-coded by rgo, the age at which 90% of the 
total stellar mass formed. Galaxies not in our sample are indicated in grey. 
Point sizes are proportional to the half-light radii. As shown by the cumula¬ 
tive distributions in the side panels, our sample provides reasonable represen¬ 
tation of all LG dwarfs in terms of proximity to a massive host and present 
day stellar mass. We are primarily missing the lowest mass systems located 
around the MW and the majority of the M31 companions, but do have some 
representation in both regimes. In terms of rgo, there are broad trends such 
as lower mass galaxies typically have large values of rgo and more distant 
galaxies tend to not be quenched. We discuss these trends in more detail in 
§3. 

ing. This was done to pick out the exact 90% value for each 
galaxy and avoid abrupt changes in the cumulative SFH that 
would compromise homogeneity in computing rgo from the 
native SFHs. For example, suppose the SFH abruptly jumps 
from 85% in one bin to 98% in a younger bin with a lull in 
star-formation in between. Without interpolation we would 
be forced to pick the value closest to 90%, whereas the inter¬ 
polation allows to consistently select a value of rgo across all 
galaxies. The assumption underlying the interpolation is that 
the SFH is constant across each bin. However, the error intro¬ 
duced by this assumption is small compared to other sources 
of uncertainty (e.g., systematics). 

We plot rgo as a function of mass and distance from the 
nearest central galaxy in Figure 1. This plot shows many 
of well-known trends: currently star-forming galaxies are lo¬ 
cated at large distances from the MW or M31 (e.g., Grcevich 
& Putman 2009), lower mass systems tend to be preferen¬ 
tially quenched compared to higher mass systems (e.g., Tol¬ 
stoy et al. 2009). However, there appears to be interesting 
scatter in both trends, particularly when the values of rgo are 
considered. We discuss these empirical relationships in more 
detail in §3. 

2.3. Sample Completeness 

Sample completeness plays an important role in assessing 
quenching in low-mass systems. Our sample is neither com¬ 
plete nor unbiased, as it was assembled based on the avail¬ 
ability of HSTAVFPC2 archival imaging. However, as shown 
in Figure 1 , the sample is reasonably representative of the LG 
dwarf population in terms of stellar mass and proximity to a 
massive host. Specifically, our sample is > 80% complete 
for systems with ^10^ Mq, but only 10-40% complete 
for galaxies with 10^ < ^10^ Mq. Similarly, our sam¬ 

ple contains ^ 50% of the galaxies known to exist within the 
virial radius of the MW or M31, and 80% of those located in 


the field. Thus, while our sample has variable completeness 
in both mass and proximity to a host, it spans the entire range 
of properties of the LG dwarf population (dashed lines on the 
side panels). In instances where completeness may compete 
with physical interpretation of the data we qualitatively dis¬ 
cuss how it may affect our results. 

3. TIME EVOLUTION OF THE QUENCHED GALAXY FRACTION 

Having defined a proxy for the quenching epoch, we now 
examine the relationship between rgo and physical properties, 
such as stellar mass and proximity to a massive host. 

3.1. Quenched Fraction as a Function of Stellar Mass and 
Lookback Time 

In Figure 2, we plot the fraction of quenched galaxies as 
a function of lookback time in four equally spaced logarith¬ 
mic bins of present day stellar mass (log(M^/M 0 )= 4 .0-5.0, 
5.0-6.0, 6.0-7.0, >7.0; hereafter we refer to these mass bins 
by their unweighted mean masses). We have defined the 
quenched fraction to be the fraction of galaxies that quenched 
at earlier an epoch than a specified lookback time (i.e., rgo> 
t). The shaded envelopes reflect the narrowest 68 % confi¬ 
dence interval of the total uncertainty (i.e., random and sys¬ 
tematic, which are inclusive of Poisson variance) that includes 
the best fit value of rgo (traced by the solid line). These un¬ 
certainty estimates are slightly conservative (i.e., too large) 
as they include possible co-variances in the SFHs that were 
derived from CMDs of similar depths, as discussed in Weisz 
et al. (2014a). 

Figure 2 shows a clear trend for lower mass dwarfs: they 
have a larger quenched fraction than more massive dwarfs at 
every epoch. By ^10 Gyr ago, ~50% of the lowest mass sys¬ 
tems (log((M^/]V[ 0 ))= 4.5) were quenched, which is at least a 
factor of 5 higher than the quenched fraction in the other three 
mass ranges considered. The quenched fraction in the lowest 
mass bin steadily grows until ~4-5 Gyr ago, by which time 
the quenched fraction has reached unity. 

In the second lowest mass bin (log((M^/M 0 ))= 5.7), the 
quenched fraction grows steadily from ^ 10% at ^ 10 Gyr 
ago to ^ 70% at ^ 3 Gyr ago. At higher masses, the 
quenched fraction is even lower. In the second highest mass 
bin (log((M^/M 0 ))= 6 . 6 ), the quenched fraction reaches 50% 
at ^ 2.5 Gyr ago. In the highest mass bin (log((M^/M 0 ))= 
8.0), the quenched fraction at the present is only ^ 40%. 

The temporal evolution of the quenched fraction can also 
provide clues as to how quickly dwarf galaxies can undergo 
quenching. For example, in the lowest mass systems, nearly 
40% appear to be quenched by ^ 12 Gyr ago, which implies 
a quenching timescale of < 1.5-2 Gyr. Such a timescale is 
consistent with quenching due to reionization or early-time 
environmental processes. Interestingly, the quenched frac¬ 
tion of the lowest mass galaxies continues to grow toward 
the present, indicating that not all low-mass systems were 
equally affected by quenching mechanisms that operated ex¬ 
clusively in the early universe (e.g., reionization). In contrast, 
galaxies in all other mass ranges have low quenched fractions 
at comparably early epochs. This broad trend is consistent 
with predictions from some reionization models, in which 
only the lowest mass systems were significantly affected by 
reionization (e.g., Ricotti & Gnedin 2005). However, other 
models suggest that the infiuence of reionization on low-mass 
galaxies is more complex than permanently quenching star- 
formation ^ 12 Gyr ago (e.g., Bemtez-Llambay et al. 2014; 
Onorbe et al. 2015). We refer the reader to Weisz et al. 
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Figure 2. The quenched fraction of LG dwarfs plotted as a function of 
present day stellar mass and lookback time. The solid lines trace the most 
likely quenched fraction values, while the uncertainty envelopes represent the 
the 68% confidence interval in the quenched fraction, including both system¬ 
atic and random uncertainties. The best fits show a broad trend in quenched 
fraction as a function of mass. The quenched fraction is higher for lower 
mass galaxies at nearly all epochs considered. Further, as a population, based 
on the quenched fraction, higher mass galaxies tend to have quenched later 
compared to lower mass galaxies. 

(2014b) for an in-depth discussion of signatures reionization 
in the SFHs of low-mass galaxies. 

In the two intermediate mass ranges, the nearly uniform in¬ 
crease in quenched fraction toward the present provides lit¬ 
tle information about the quenching timescale. In contrast, 
the highest mass galaxies are more informative for about the 
rapidity of quenching. In this mass range (log((M^/MQ))= 
8.0), the quenched fraction grew from 0 to 40% between 2 
and 4 Gyr ago. This rapid rise indicates that galaxies of this 
mass can quench on fairly rapid timescales, although without 
a knowledge of the infall times we cannot determine the exact 
quenching timescale (see §4). Further, the six systems in this 
mass range may not be broadly representative of galaxies in 
this mass range, given that four are dEs in M31 group (M32, 
NGC 147, NGC 185, NGC 205) and the other two are MW 
companions Fornax and Sagittarius. The nearly synchronized 
quenching of M31 dEs may be indicative of a global event in 
the M31 group (e.g., simultaneous satellite accretion, major 
merger; Fardal et al. 2008; Hammer et al. 2010), which may 
detract from the generality of this finding. For a further dis¬ 
cussion of quenching in M31 satellites, we refer the reader to 
§3.3 and Weisz et al. (2014c). 

Up to this point, we have not considered the role of uncer¬ 
tainties and selection effects in our results and interpretation. 
The uncertainties on the quenched fraction are most notice¬ 
able for the lowest mass galaxies, which typically show broad 
68% confidence intervals on their values of rqq. This breadth 
is the results of a small number of stars on the CMDs of the 
lowest mass systems, as well as the presence of BSs, which 
introduce a large dispersion into rqq. For the more massive 
systems, which have more populated CMDs, the effects of 
uncertainties in the quenched fraction are noticeably less pro¬ 
nounced. 

Completeness of the sample can also affect our findings. 
Many of the faint MW galaxies that are not in our sample 


appear to have quenched at lookback times of > 10-12 Gyr 
ago (e.g., de Jong et al. 2008; Sand et al. 2009, 2010, 2012; 
Brown et al. 2012, 2014; Okamoto et al. 2012). A naive 
completeness correction, in which we assume that all miss¬ 
ing galaxies are similar to the faint MW companions, implies 
that the quenched fraction could be as high as ^80% by 10-12 
Gyr ago. However, as demonstrated in Weisz et al. (2014c), 
not all quenched low-mass systems in the LG stopped form¬ 
ing stars >10 Gyr ago. Deep imaging of Andromeda XVI 
(M^ ^10^ Mq) reveals a quenching epoch of ^ 5 Gyr ago. 
Similarly, other low-mass galaxies in our sample such as An¬ 
dromeda XH and Leo T have significant amounts of star for¬ 
mation <10 Gyr ago, even when full uncertainties are con¬ 
sidered. Thus, on the whole, the simple assumption that all of 
the lowest-mass quenched galaxies are similar to the faintest 
MW satellites may not be an accurate extrapolation. Finally, 
we note our sample is 80% complete for M^ > 10® M©, indi- 
eating that incompleteness in this regime is not as important 
as for lower mass systems. 

3.2. Quenched Fraction as a Function of Environment 

In Figure 3, we plot the quenched fraction as a function 
of lookback time in bins of current distance from a massive 
host. We have divided the sample into four equally spaced 
radial bins in multiplies of Rviriab ranging from R<0.5 Ryiriai 
to R> 1.5 RyiriaL where we adopt Ryiriai=300 kpc for both M31 
and the MW (McConnachie 2012). 

From this plot, it is clear that environment plays an im¬ 
portant role in quenching low-mass galaxies. At the present 
day, few dwarfs located outside Rymai are quenched, while 
those located inside Rymai are almost entirely quenched. This 
trend is essentially a re-statement of the long-established 
LG morphology-density relationship (e.g., Mateo 1998; van 
den Bergh 2000; Grcevich & Putman 2009, and references 
therein). 

For dwarfs located inside Ryiriai, the temporal evolution of 
the quenched fraction is essentially independent of the cur¬ 
rent proximity to M31 or the MW, and is statistically indis¬ 
tinguishable in both bins within R<300 kpc. The small dif¬ 
ferences at recent times are due to the presence of gas-rich 
galaxies LGS 3 and IC 10 inside the virial radius of M3L 
The independence between quenching time and current prox¬ 
imity to a massive host is consistent with a picture in which 
quenching occurs when or soon after a galaxy enters the virial 
radius. This conclusion does not depend on selection effects 
and uncertainties in the values of rqq. 

In principle, the SFHs may also provide clues to the planar 
satellite structures that are known to exist around the MW and 
M31 (e.g., Lynden-Bell 1976; Kroupa et al. 2005; Ibata et al. 
2013). For example, if the planes were accreted as a group 
at the same time, the process of accretion may leave system¬ 
atic signatures in the quenching epochs that are different from 
off-plane systems. Unfortunately, we do not detect any such 
clear trends. In the case of M31, only two of our galaxies. 
And H and And V are off-plane, limiting any general com¬ 
parison. Similarly for the MW sample, there is no clear dif¬ 
ference in the SFHs of planar and non-planar satellites. The 
lack of a correlation does not rule out any large scale accre¬ 
tion scenario, as dynamical mixing may have erased the initial 
conditions of any such event. 

3.3. Quenching Properties of the Milky Way and MSI 

Satellites 
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Figure 3. The same as Figure 2, except for proximity to M31 or the MW 
instead of stellar mass. Galaxies inside the virial radius of M31 or the MW are 
largely quenched (the exceptions are IC 10 and LGS 3), while those outside 
the virial radius are not (akin to the LG morphology-density relationship). 
However, within 300 kpc there is no clear trend between present day distance 
from a host and quenching epoch. 

As discussed in §3.1, there are hints that the M31 and MW 
satellites may have different quenching characteristics. We 
compare these two systems in detail in Figure 4. Here, the 
quenched fractions show qualitatively different trends. While 
the quenched fraction for the MW satellites steadily increases 
toward the present, the quenched fraction for M31 dSphs 
jumps from 20% at ^ 7 Gyr ago to 100% by ^ 5 Gyr ago. 
Further, the M31 dEs appear to have rapidly quenched be¬ 
tween ^ 2 and 4 Gyr ago, a different timescale than both sets 
of dSphs. 

Taken at face value, these results suggest that satellites in 
the two sub-groups may have systematically different quench¬ 
ing characteristics, despite the MW and M31 having similar 
stellar masses. The broader implication is that satellite be¬ 
havior may be sensitive to the specific accretion history of 
the host galaxy, as opposed to simply scaling with host mass. 
However, our data only provide hints of these differences due 
to large uncertainties in the SFHs of M31 satellites and in¬ 
completeness in both the MW and M31 satellite samples. 
Indeed the uncertainties allow there to be identical quench¬ 
ing timescales between the two systems. Upcoming HST 
programs aimed at systematically imaging the M31 satellites 
(GO-13739, PLSkillman; GO-13699, PI: Martin) should pro¬ 
vide the data necessary to make more definitive statements 
about the quenching epochs of M31 satellites (e.g. Weisz et al. 
2014c). 

3.4. The Quenched Fraction over 7 dex in Stellar Mass 

The fraction of quenched galaxies as a function of redshift 
and stellar mass is widely used as an observational constraint 
for modeling quenching processes in groups and clusters of 
galaxies (e.g., Mandelbaum et al. 2006; van den Bosch et al. 
2008; Peng et al. 2010; Bauer et al. 2013; Moustakas et al. 
2013; Wetzel et al. 2013). However, the majority of such 
studies have been limited to galaxies with M^ > 10^-10^^ 
Mq, since fainter galaxies are extremely challenging to de¬ 
tect beyond the very nearby universe. As a result, there is no 


Redshift (z) 

10 5 2 1 0.5 0.2 0.1 0 



Lookback Time (Gyr ago) 

Figure 4. Time evolution of the quenched fraction of MW dSphs (red), M31 
dSphs (blue), and M31 dEs (dashed line; shown only for qualitative refer¬ 
ence). The most likely quenched fractions (points) appear different for the 
M31 and MW dSphs. The quenched fraction of M31 dSphs increased rapidly 
at intermediate ages (5-7 Gyr ago), while the quenched fraction of MW dSphs 
has steadily grown over the history of the universe. These contrasting pattens 
may provide clues to differences in the evolution of the sub-groups, however 
the uncertainties are too large to draw definitive conclusions. In comparison 
to the dSphs, the quenched fraction of the M31 dEs increased rapidly between 
2 and 3 Gyr ago. 

unified picture of quenching properties and processes over a 
full compliment of galaxy masses and lookback times. 

At the same time, the observational evidence that quench¬ 
ing in lower mass galaxies is strongly driven by environment 
is steadily growing. For example, Geha et al. (2012) show the 
fraction of quenched galaxies with M^ ^ 10^-^ Mq are rare at 
distances >1.5 Mpc from a massive host. Within the near¬ 
est few Mpc, KKR 25 and KKs3 are only known examples 
of quenched, isolated dwarf galaxy (e.g, Weisz et al. 2011; 
Makarov et al. 2012; Karachentsev et al. 2015). The pre¬ 
sumed^ ^ The rarity of quenched, low-mass field galaxies sug¬ 
gests that most low-mass galaxies quenched in higher-density 
environments, although it is not possible to rule out complete¬ 
ness effects in the current surveys. 

Therefore, exploring a broader picture of quenching re¬ 
quires comparing galaxies with respect to both mass and en¬ 
vironment. In the remainder of the section, we compared 
quenching in the LG dwarfs with higher mass galaxies from 
all environments. In the next section, we discuss the LG 
dwarfs alongside satellite populations to assess the specific 
role of environment. 

In Figure 5 and Table 2, we illustrate the evolution of the 
quenched fraction from z = 0 -1 (^ 0-8 Gyr ago) for galaxies 
spanning a range of 7 dex in stellar mass. In the top panel, we 
have combined the quenched fractions derived in this paper 
for galaxies with M^ < 10^ Mq from Wetzel et al. (2013) (M^ 
^ 10^^ < 10^^^ Mq; z = 0 -1). This figure provides a 

first view of the quenched fraction over large ranges in stellar 
mass and redshift. Note that for consistency with the analy- 

It is not possible to rule out completeness effects in the current low-mass 
galaxy surveys. However, any attempts to quantify the numbers and charac¬ 
teristics of undiscovered low-mass field galaxies would be pure speculation. 
We therefore take the current results at face value while acknowledging the 
competing effects of completeness. 
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Figure 5. The evolution of the quenched fraction of all galaxies over 
~ 7 dex in stellar mass from z = 0 - 1. Top-The quenched fraction 
for low mass galaxies (M* <10^ M©) are from this paper, and the 
higher mass systems (M* > 10^ M©) are from Geha et al. (2012) and 
Wetzel et al. (2013). bottom - The thin colored lines are the same 
as in the top panel, while the thick lines are the quenched fraction of 
all and field galaxies from Geha et al. (2012) and Karachentsev et al. 
(2013), as described in §3.4. The envelopes for the Karachentsev 
et al. (2013) data shows the range of quenched fractions for whether 
dTrans are considered quenched or not. In both panels, the “U” shape 
in the the quenched fraction as a function of mass may refiect the 
mass scales at which different quenching mechanism operate. At 
low masses, external processes (e.g., reionization, ram pressure, tidal 
stripping) are thought to drive quenching, while quenching in higher 
mass systems is controlled by halo mass (e.g., gas accretion rates, 
AGN feedback, virial shocks, stellar feedback). The ubiquitously 
low quenched fraction at intermediate masses (M* > 10^-10^^ M©) 
may indicate a transition region in which neither type of quenching 
mechanism is particularly effective. 


sis of Wetzel et al. (2013), we compute the LG dwarf stellar 
masses at the specified redshift. There are generally negligi¬ 
bly different than the z = 0 masses used elsewhere in the paper, 
as most dwarfs only formed a small fraction of their mass at 
z<l. 

The most notable feature in this plot is the minimum of 
the quenched fraction as a function of mass, at a given red- 
shift. At each redshift considered, the quenched fraction 


Table 2 

Quenched Fraction of Galaxies as a Function of Redshift and 
Mass 


log(M*(z)) 

(Mo) 

(1) 

Quenched 

Fraction 

(2) 

log(M*(z)) 

[(Mo) 

(3) (4) 

Quenched 

Fraction 

(5) 

z = 0 


z = 0.2 


4.5 

1.00 

4.5 

1.0 

5.7 

0.73 

5.7 

0.73 

6.6 

0.56 

6.5 

0.44 

8.0 

0.50 

7.9 

0.25 

II 

o 


z = 0.6 


4.5 

0.83 

4.5 

0.67 

5.7 

0.64 

5.5 

0.55 

6.5 

0.33 

6.5 

0.2 

7.9 

0.00 

7.9 

0.00 

II 

p 

bo 


z=1.0 


4.5 

0.57 

4.5 

0.57 

5.6 

0.40 

5.6 

0.3 

6.4 

0.2 

6.4 

0.1 

7.8 

0.00 

7.8 

0.00 


Note. — The fraction of quenched LG dwarf galaxies as a 
function of stellar mass and redshift. This table corresponds 
to the values plotted in Figures 5 and 6. We emphasize that 
due to selection effects, these values should be considered in a 
qualitative sense, as discussed in §3.4. 

Steadily declines from maxima at the highest and lowest stel¬ 
lar masses toward a minimum at 10^-10^^ M 0 . The “U" 
shape in the fraction of quenched galaxies may provide clues 
for which quenching mechanisms operate at different mass 
scales, which we now discuss in more detail. 

At the highest masses (M^>10^^ Mq), quenching may be 
driven by several internal processes including AGN feedback, 
virial shock heating, restricted gas accretion, stellar feedback 
and/or disk instabilities (e.g., Silk & Rees 1998; Keres et al. 
2005; Birnboim et al. 2007; Dekel et al. 2009; Woo et al. 
2013). Although the relative contribution of each mechanism 
to quenching remains unclear, the broad picture is that the 
depth of the potential well (i.e., halo mass) largely drives the 
initial quenching by depleting gas and suppressing cooling. 
Further, star formation may be prevented by a combination 
of mechanisms (such as those above) or via other pathways 
such as heating from asymptotic giant branch stars (e.g., Con¬ 
roy et al. 2014) that appear incapable of initial quenching, but 
may provide sufficient heating to maintain a massive galaxy’s 
quenched state. 

At slightly lower masses (M^^IO^^ Mq), halo mass may 
play a diminished role in quenching. The ratio of stellar-to- 
halo mass reaches a global maximum at M^~ 10 ^^-^ MQ(e.g., 
Moster et al. 2010; Behroozi et al. 2013), suggesting that stel¬ 
lar feedback may play a more important role in quenching in 
this mass range (e.g., Moster et al. 2010; Hopkins et al. 2011, 
2012, 2014). 

The next lowest mass data points (M^<10^ Mq) are from 
the LG sample. Compared to the Wetzel et al. (2013) data, 
these galaxies are (a) highly dark matter dominated and (b) 
represent an obvious transition to a group environment. As a 
result, qualitatively different quenching mechanisms may be 
at work. Typically, quenching of galaxies in this mass range is 
ascribed to environmental processes such as ram pressure or 
tidal stripping (e.g., Mayer et al. 2001, 2006), both of which 
require the presence of a massive galaxy, i.e., a group envi¬ 
ronment, to be effective. Putative environmentally indepen¬ 
dent quenching mechanisms such as stellar feedback and gas- 
consumption due to star formation (e.g., Dekel & Silk 1986), 
are not highly viable quenching mechanisms. Stellar feed- 
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Table 3 

Quenched Fraction For the Nearby Galaxy Catalog Karachentsev et al. (2013) 


(log(M./M0)> 

Quenched Fraction 

Nrotal 

Nquenched 

dTrans= 

Quenched Fraction 

Nquenched 
dTrans = 

(1) 

(2) 

(3) 

star-forming 

(4) 

(5) 

not star-forming 
(6) 

All Galaxies 

5.68 

0.60 

15 

9 

0.67 

10 

6.74 

0.35 

139 

49 

0.50 

70 

7.60 

0.13 

310 

38 

0.15 

47 

8.55 

0.04 

201 

8 

0.04 

8 

9.53 

0.04 

94 

4 

0.04 

4 

10.32 

0.19 

31 

6 

0.19 

6 

Field Galaxies 

5.80 

0.67 

6 

4 

0.67 

4 

6.77 

0.19 

74 

14 

0.33 

25 

7.60 

0.05 

237 

12 

0.07 

16 

8.55 

0.03 

164 

5 

0.03 

5 

9.51 

0.01 

80 

1 

0.01 

1 

Satellite Galaxies 

5.58 

0.56 

9 

5 

0.67 

6 

6.71 

0.54 

65 

35 

0.69 

45 

7.60 

0.36 

72 

26 

0.43 

31 

8.53 

0.08 

37 

3 

0.08 

3 

9.60 

0.21 

14 

3 

0.01 

3 


Note. — The z = 0 quenched fraction of galaxies as determined from analysis of the Karachentsev et al. (2013) 
nearby galaxies catalog (D < 11 Mpc). The quenched fractions were computed by comparison of listed mor¬ 
phological types and 3D distances as described in §3.4. Due to potential confusion between dSph and dTrans 
morphological types at low stellar masses, we have included quenched fractions in which dTrans are counted star¬ 
forming (columns (4) and (5)) or are counted as quenched (columns (5) and (6)). LG galaxies have been exulted 
from this analysis. 


back creates outflows of hot gas, but is not energetic enough to 
completely unbind cold gas from dark matter halos (e.g., Mac 
Low & Ferrara 1999; Govemato et al. 2010, 2012). Similarly, 
gas consumption from star formation is not a likely quenching 
mechanism. Gas consumption timescales for most gas-rich 
dwarfs are tens to hundreds of Gyr (e.g., van Zee et al. 1997; 
Skillman et al. 2003), and, if acting alone, would produce few 
quenched dwarfs in a Hubble time. 

At the lowest masses (M^<10^ Mq), reionization may po¬ 
tentially playing a role in quenching (e.g. Efstathiou 1992; 
Bullock et al. 2000; Ricotti et al. 2002). Although often 
treated as a global, and thus an environmentally independent 
process, detailed simulations suggest that reionization is not 
necessarily uniform in time, intensity, and/or environment; 
therefore its impact on low mass systems may be highly vari¬ 
able (e.g., Busha et al. 2010; Ocvirk et al. 2013). As discussed 
in Weisz et al. (2014b) and in §3.1, large variability in the 
SFHs of the lowest mass galaxies may be consistent with an 
inhomogenous reionization scenario. 

Unfortunately, data from the LG and Wetzel et al. (2013) 
provide no sampling in the minimum mass range (M^^ 10^- 
10^^ Mq). However, nearby galaxy surveys can All in this 
range at z = 0 providing some sampling for the behavior of 
these systems. In the bottom panel of Figure 5 we add the 
quenched fractions for all and held galaxies from the SDSS- 
based study of Geha et al. (2012) and from the nearby galaxy 
catalog of Karachentsev et al. (2013). For the former dataset, 
held galaxies are those with projected distances >0.5 Mpc 
from a massive host. 

For the latter, we have estimated the quenched fraction from 
the nearby galaxy catalog (1<D<11 Mpc) of Karachentsev 
et al. (2013) as follows. We first convert all listed M^ values 
into stellar masses assuming MQfLQ=l and a Mb^q =5.45. 
We designate all galaxies with M^ >10^^ Mq to be ‘central’ 
galaxies, and compute the 3D distances (based on the listed 
line of sight distances and sky coordinates) between each cen¬ 


tral and the other 861 galaxies in the catalog. Galaxies that are 
within 0.5 Mpc of a central are considered to be satellites, and 
more distant non-centrals are designated as held galaxies. We 
then compute the quenched fraction by comparing listed mor¬ 
phologies; those with T<0 are considered to be quenched and 
those with T>0 are not. At M^ < 10^ Mq, the Karachentsev 
et al. (2013) suffers from decreasing completeness, and the 
authors also note that there is some ambiguity in morpholog¬ 
ical typing between star-forming transition objects that have 
faint GALEX fluxes but no detectable Ha (i.e., dTrans) and 
gas-free, non-star forming dSphs. To mitigate this latter ef¬ 
fect, we have computed the morphological-based quenched 
fractions in both cases (i.e., all are dTrans are quenched or 
all dTrans are star-forming). These ranges are reflected by 
the envelops in the bottom panel of Figure 5. Throughout 
this process, all LG galaxies were excluded from considera¬ 
tion. As discussed below, the quenched fractions determined 
by this method may not be as quantitative as that of Geha et al. 
(2012), but this does provide an extension to lower masses 
that is not possible with SDSS. As a sanity check, the mor¬ 
phological approach provides very good overall agreement 
with expectations from both LG and ANGST SFHs. We list 
the quenched fractions derived from the Karachentsev et al. 
(2013) in Table 3. 

The addition of literature data at z = 0 nicely Alls in missing 
mass range from the top panel of Figure 5 and helps paint a 
picture of quenching trends over a wide range of stellar mass. 
For the intermediate mass range (M^ ~ 10^ - 10'° M 0 ), the 
quenched fractions from Geha et al. (2012) and Karachentsev 
et al. (2013) are nearly zero. Both exhibit a small to modest 
upturn at M^ > 10^^ Mq, which is qualitatively consistent 
with the Wetzel et al. (2013) results. 

At M^ < 10^ Mq, the Karachentsev et al. (2013) quenched 
fractions show a signiflcant upturn. This is particularly in¬ 
teresting for the held sample, as is it not clear what mecha¬ 
nisms can effectively quench such galaxies in the absence of 
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Figure 6. The evolution of the quenched fraction of satellite galaxies over ~ 
7 dex in stellar mass. The LG data is the same as in Figure 5, while the liter¬ 
ature values represent non-LG galaxies by Geha et al. (2012); Karachentsev 
et al. (2013); Wetzel et al. (2013); Slater & Bell (2014); Wheeler et al. (2014); 
Phillips et al. (2015). The grey envelope for the Karachentsev et al. (2013) 
data shows the range of quenched fractions for whether dTrans are considered 
quenched or not. For M* <10^ Mq, there is a clear increase in quenched 
fraction toward lower masses, whereas above this mass the quenched fraction 
is roughly constant, although the analysis of Slater & Bell (2014) suggesting 
an upturn is possible. This suggests that satellites with M* >10^ Mq are not 
resistant to quenching, while those of decreasing mass are more vulnerable to 
quenching mechanisms that operate in high density environments (e.g., ram 
pressure, tidal stripping) and possibly reionization at the lowest masses. 

environment. One possibility is that there are so-called ‘back- 
splash’ galaxies, which had an interaction with a massive host 
and have since been ejected outside of the virial radius, plac¬ 
ing them into the quenched field galaxy category. Similarly, 
analysis of the ELVIS Nbody simulations (Garrison-Kimmel 
et al. 2013) by Deason et al. (2014) suggest that dwarfs out¬ 
side the virial radius of a massive galaxy have a moderately 
high merger rate, which may be capable of inducing quench¬ 
ing. Alternatively, it may be the case that lower-mass sys¬ 
tems are more vulnerable to gas-loss due to feedback than 
current simulations suggest, which would introduce a mass- 
dependent quenching trend in isolated systems. Suggestions 
that dwarfs do not follow the stellar-halo mass relation from 
higher mass galaxies (e.g., Garrison-Kimmel et al. 2014), i.e., 
they more form more stars for a given halo mass, may also 
play a role in this possibility. Taken at face value, these data 
indicate that there are a modest number of quenched low-mass 
galaxies in the nearby universe, although there does not re¬ 
main a clear physical reason for the observed quenched frac¬ 
tion in the field particularly as the majority of galaxies in this 
sample are too massive to have been completely quenched by 
photo-heating due to cosmic reionization (e.g., Bullock et al. 
2000; Ricotti & Gnedin 2005). 

Aside from physical effects, we must also consider the pos¬ 
sibility of systematic uncertainties in the Karachentsev et al. 
(2013) data. In particular, inaccurate morphological classi¬ 
fications at lower masses, incompleteness, and/or inaccurate 
distances could all infiuence the quenched fraction beyond the 
nominal Poisson errors. We address this confusion by com¬ 
puting the quenched fractions in the cases that all dTrans and 
quenched/not quenched, and plot the range as envelopes in 
Figure 5. While these results are intriguing, it is clear that 
a more detailed study of low-mass field galaxies in the very 


nearby universe will help clarify the origin of upturn in the 
Karachentsev et al. (2013) data. 

The redshift evolution in Figure 5 is consistent in a pic¬ 
ture of essentially monotonic quenching. That is, once a 
galaxy becomes quenched, it is unlikely to go back to be¬ 
ing star-forming. There are many mechanisms that can in¬ 
duce and maintain a galaxy’s quenched status (e.g., feedback, 
ram pressure), but essentially only gas accretion can re-ignite 
star formation. Given that this single process must compete 
with a plethora of quenching mechanisms, it is unlikely to 
observe significant decreases in the quenched fraction toward 
the present. Indeed, this trend is observed in both the LG and 
Wetzel et al. (2013) data. 

Applying this same reasoning to the z = 0 literature data, we 
can conclude that galaxies with 10^-10^^ M© have had 
similarly low quenched fractions in the past. That is, galaxies 
in this mass range truly appear to be resistant to quenching, 
irrespective of their environment. 

3.5. Satellite Quenching 

Although we have just discussed the LG in the context of all 
galaxy quenching (i.e., centrals and satellites), a more conven¬ 
tional approach is to analyze only centrals or only satellites. 
Thus, in Figure 6, we plot the redshift evolution of our LG 
data along with satellite quenched fractions from several z = 0 
satellite quenching studies. 

The overall finding of these studies is that satellites with 
^10^ M© have nearly constant quenched fractions over 
over a large range of stellar mass (M^ ^ 10^ - 10^^ M©). 
There are some minor inconsistencies in trends reported by 
Slater & Bell (2014) and Wheeler et al. (2014), which ap¬ 
plied different selection criteria to the same Geha et al. (2012) 
dataset. This disagreement highlights another challenge in 
quantitatively interpreting data on quenching from different 
studies, a point that can be extended to our derivation of the 
satellite quenched fraction from the Karachentsev et al. (2013) 
catalog. 

Despite subtle catalog differences, the general belief is that 
satellites with ^10^ M© are fairly resistant to quenching, 
even when acted on by ram pressure and/or tidal effects. The 
exact timescales of quenching are still sensitive to choices in 
running and interpreting simulations. However, commonly 
referenced timescales for becoming quenched after a satel¬ 
lites enters the virial radius of a massive galaxies are >5-9 
Gyr (Slater & Bell 2014; Wheeler et al. 2014; Phillips et al. 
2015). Although Slater & Bell (2014) suggest that follow¬ 
ing a pericentric passage, quenching must have more rapid 
timescales of ^ 1-2 Gyr, which is consistent with our current 
understanding of Leo I. We discuss more details on the rela¬ 
tionship model-based infall times and SFH-based quenching 
epochs in the next section. 

For galaxies with ^10^ M©, the increase in quenched 
fraction with decreasing mass suggests that lower mass sys¬ 
tems are more vulnerable to the environment quenching ef¬ 
fects in a group environment. However, a comparison with the 
Karachentsev et al. (2013) field data points in Figure 5 sug¬ 
gests that secular processes (e.g., feedback) may also be im¬ 
portant. It is unclear how to reconcile these two figures in light 
of the current paradigms about gas-loss in low-mass galax¬ 
ies, which do not make a clear prediction for mass-dependent 
quenching in the field. While this observation may hint at new 
insight into the physics of quenching in low-mass galaxies, it 
also is clouded by potential selection effects and small num¬ 
ber statistics as discussed in §3.4. It is further complicated 
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Figure 7. Constraints on the delay time (the difference between infall time 
and the quenching epoch) for select galaxies associated with the MW. The 
infall times on the y-axies are based on analysis of Via Lactea presented in 
Rocha et al. (2012) and our values of rgo are on the x-axis. Point sizes are 
proportional to the half-light radii, and error bars represent the total (grey) 
and random uncertainties (navy; in the case of rgo). Galaxies located to the 
upper left quenched prior to infall, while those in the lower right quenched 
after infalL Galaxies that quenched after infall appear to have done so within 
~ 1-4 Gyr after entering the virial radius of the MW. 


by our primitive understand of the effects of reionization on 
low-mass galaxies (e.g,. Bemtez-Llambay et al. 2014; Onorbe 
et al. 2015). While many in this sample are above the stel¬ 
lar mass range assumed to be affected by reionization, a new 
generation of low-mass galaxy simulations suggests that the 
effects of reionization may extended to higher mass systems 
(e.g., Bemtez-Llambay et al. 2014). In such cases, it may be 
that reionization in tandem with other effects (e.g., feedback) 
can quench or strongly suppress star formation in galaxies 
more massive than posited by the first generation of models 
(e.g., Ricotti & Gnedin 2005). For deeper empirical and the¬ 
oretical understanding, more extensive surveys of low-mass 
galaxies in isolated environments are needed to understand 
the significance of these trends revealed in this paper. 

4. QUENCHING VS. INFALL TIMESCALES 

Constraints on the processes responsible for quenching and 
its efficiency can come from comparing the infall time and 
quenching time. The ‘delay time’ between the two is a proxy 
for how long a galaxy can retain its gas and continue forming 
stars after entering the virial radius of a massive host. Alterna¬ 
tively, if quenching happens prior to infall, this may indicate 
that non-environmentally driven processes (e.g., reionization) 
were responsible. 

While quenching epochs come directly from our SFH mea¬ 
surements, infall times cannot be measured directly and in¬ 
stead must be inferred from models. A common approach 
is to match sub-halos from N-body simulations with observed 
galaxies by comparing their masses (via abundance matching) 
and available phase space information (e.g., tangential and ra¬ 
dial velocities). Rocha et al. (2012) provide a clear example of 
this type of modeling applied to galaxies associated with the 
MW. Using Via Lactea II (Diemand et al. 2008), they quantify 
the correlation between the present day orbital energy and the 
epoch at which a galaxy most recently entered the virial radius 


of the MW. Combined with available radial and proper mo¬ 
tions (when available), they predict the infall time probability 
distribution functions for a dozen dwarf galaxies associated 
with the MW. In the remainder of this section, we compare 
their infall time predictions with our quenching timescales 
and explore implications for quenching mechanisms. 

Figure 7 shows several interesting trends in the comparison 
between rgo and the predicated infall times. First, several of 
the lowest mass galaxies appear to have quenched before in¬ 
fall. For example, Hercules quenched ^12 Gyr ago, but only 
entered the virial radius of the MW ^ 5 Gyr ago, suggesting 
that it spent the majority of its life away from the environmen¬ 
tal influence of the MW. Six other galaxies (Leo IV, CVn I, 
CVn II, Leo II, Sculptor, and Draco) also have earlier quench¬ 
ing than infall times. Of these six galaxies, more than half 
(Hercules, Leo IV, CVn I, and CVn H) have been previously 
identified as putative ‘fossils of reionization’, i.e., their star 
formation was quenched due by reionization (e.g., Bovill & 
Ricotti 2011; Brown et al. 2012, 2014). If they were located 
outside the virial radius of the MW prior to quenching, the 
likelihood of quenching due to environmental processes (e.g., 
ram pressure) seems drastically lower, increasing the likeli¬ 
hood that reionization played a significant role in quenching 
their star formation. 

Second, at the other extreme, Carina and Fornax appear to 
have quenched several Gyr after their infall times. Figure 
7 suggests a delay of ^ 5 Gyr between when these galax¬ 
ies fell in, and when their star formation stopped. This 
scenario would be qualitatively similar to the ‘delayed-then- 
rapid’ quenching scenario for more massive satellites (M^ > 
10^ - 10^^ Mq) that was proposed by Wetzel et al. (2013). 
In this scenario a satellite continues to form stars for ~ 2-4 
Gyr after infall, before rapidly quenching in less than 1 Gyr. 
We will explore more detailed modeling of this scenario in a 
future paper. 

Third, there are a handful of galaxies that quenched near 
or shortly after their infall time. One example is Leo 1. Leo 
I appears to have fallen into the MW’s halo ^2.3 Gyr ago 
(Rocha et al. 2012; Sohn et al. 2013). ^ 1 Gyr ago, Leo I 
appears to have reached pericenter at ^ 90 kpc from the MW, 
which is coincident with its value of rgo. This suggests an 
extremely small delay time, which is in stark contrast to the 
comparably luminous MW dSph Fornax. However, Leo I is 
believed to be on a highly inclined orbit, while Fornax is more 
circular, suggesting that exact orbital trajectory may play a 
role in a the delay time between infall and quenching. 

Carina, which has a lower stellar mass than Fornax, also 
quenched ~ 5 Gyr after infall, supporting the idea that mass 
may not be the sole determinant in quenching timescale, and 
exact orbital history may matter. Interestingly, in the case 
of Carina, there is an open question as to whether its SFH 
shows multiple distinct bursts followed by periods of little or 
no star formation (e.g., Smecker-Hane et al. 1994; Hurley- 
Keller et al. 1998; Bono et al. 2010). Interpreted in the context 
of infall, it could be that Carina has made more than one pas¬ 
sage around the MW, but only lost all of its gas and became 
completely quenched within the last few Gyr. As shown in 
Weisz et al. (2014a), our SFH of Carina has similar bursts and 
lulls in its SFH to previous literature derivations. However, 
when all uncertainties are considered, it is possible that the 
SFH was never completely quiescent between bursts. 

Finally, Rocha et al. (2012) suggest that Leo T only en¬ 
tered MW virial radius within the last Gyr. Observationally, 
the impact of this recent accretion appears to be minimal. 
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Leo T has a large baryonic gas fraction (0.8) and displays spa¬ 
tially symmetric HI and stellar distributions (e.g., Irwin et al. 
2007; Ryan-Weber et al. 2008). Further, Leo T has nearly 
constant star formation over the last few Gyr, showing little 
evidence for Leo I-like enhancement upon infall. It does lack 
detectable Ha emission, which could be interpreted as a break 
in the SFH. However, it could also be due to stochastic sam¬ 
pling of the stellar IMF and/or cluster mass function, which 
is believed to be common in such low-mass systems (e.g., Fu- 
magalli et al. 2011; da Silva et al. 2012). 

5. SUMMARY 

We have used the SFHs of 38 LG dwarf galaxies to explore 
the quenching characteristics of low-mass galaxies (10^ < 
<10^ Mq). The SFHs were uniformly measured from the 
analysis of resolved star CMDs based on HST/WFPC2 imag¬ 
ing that were presented in Weisz et al. (2014a). Using a proxy 
for the quenching epoch as the lookback time which 90% of a 
galaxy’s stellar mass formed, we found the following results: 

1. Lower mass galaxies tend to quench prior to higher 
mass galaxies. Specifically, the quenched fraction of 
the lowest mass galaxies (log((M^/M0))=4.5) is 0.4 
at ^ 10 Gyr ago and 1.0 at ^ 3-4 Gyr ago. In con¬ 
trast, at the same epochs, the highest mass systems 
(log((M*/MQ))= 8.0) have a quenched fraction of < 
10% and ^ 40%. The quenched fraction appears to 
vary smoothly as a function of mass and lookback time 
between these extremes. Due to incompleteness of our 
sample, the quenched fraction of the lowest mass galax¬ 
ies may be larger at early times. However, several of 
the lower mass dSphs in the LG show evidence for 
extended SFHs, so the assumption that all low-mass 
galaxies quenched early may not be entirely accurate. 

2. Within Rviriai of either the MW or M31, there is lit¬ 
tle correlation between present distance to a host and 
quenching epoch. Outside of Ryiriai the quenched frac¬ 
tion is effectively zero at all lookback times. 

3. There are hints of systematic differences in the quench¬ 
ing epochs of M31 and MW dSphs. The quenched frac¬ 
tion of M31 dSphs jumps from 0.2 at ^ 7 Gyr ago to 1.0 
by r\j 5 Gyr ago, whereas the quenched fraction of MW 
dSphs steadily increases over cosmic time. Further, the 
four M31 dEs all appear to have quenched between ^ 
3 and 4 Gyr ago. However, the sample size and shallow 
CMDs of the M31 companions make it challenging to 
draw any definitive conclusions. 

4. There are clear trends in the quenched fraction over ^ 
7 dex in stellar mass. The lowest (M^ ^10^ Mq) and 
highest mass (M^ ^ 10^^-^ Mq) galaxies always have 
the highest quenched fractions between z = 0 -1, while 
galaxies with M^ ^10^-10^^ Mq always have the low¬ 
est quenched fraction. Given that the quenched fraction 
essentially monotonically increases toward the present, 
we conjecture that galaxies with M^ ^ 10^ - 10^^ Mq 
are universally the most difficult to quench. 

5. A comparison with the predicted infall times from 
Rocha et al. (2012) for select galaxies associated with 
the MW allows us to estimate delay times (i.e., the time 
between infall and quenching). We find that lower mass 
systems tended to have quenched prior to infall, while 


higher mass systems quenched after infall, and exhibit 
delay times of rsj - 4 Gyr. 
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